Quenched binary Bose-Einstein condensates: spin domain formation and coarsening 
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We explore the time evolution of quasi- ID two component Bose-Einstein condensates (BEC's) 
following a quench from one component BEC's with a U(l) order parameter into two component 
condensates with a U(l) x Z2 order parameter. In our case, these two spin components have a 
propensity to phase separate, i.e., they are immiscible. Remarkably, these spin degrees of freedom 
can equivalently be described as a single component attractive BEC. A spatially uniform mixture 
of these spins is dynamically unstable, rapidly amplifing any quantum or pre-existing classical spin 
fluctuations. This coherent growth process drives the formation of numerous spin polarized domains, 
which are far from the system's ground state. At much longer times these domains grow in size, 
coarsening, as the system approaches equilibrium. The experimentally observed time evolution is 
fully consistent with our stochastic-projected Gross-Pitaevskii calculation. 

PACS numbers: 75.75.+a,75.40.Gb 



Ultracold atomic gases are unique systems for studying 
phase transitions where the full range from adiabatic to 
diabatic can be easily accessed in the laboratory. A prime 
example of this is the transition from superfluid (SF) to 
Mott-insulator (MI) in an optical lattice: when the tran- 
sition is crossed slowly, a nearly T = SF transforms 
into a nearly T = MI [1 j; however, when the system is 
quenched by rapidly entering the MI regime, it exhibits 
rapid dynamics before dephasing into a highly excited, 
high temperature final state [2H3]. Here we study a sim- 
ilar quantum quench in a two component spinor BEC, 
where the spin degree of freedom is initialized in a highly 
excited state. We follow the resulting rapid dynamics 
during which spin domains first form, and subsequently 
slowly relax towards equilibrium as the domain size in- 
creases and the domain number decreases (see Fig [I]). 
To our knowledge, this is the first observation of domain 
coarsening py in a spinor BEC. 

The establishment of out of equilibrium domains 
formed by quenching through a phase transition is ubiq- 
uitous in physical systems ranging from grain formation 
in minerals [6] , structure growth in the early universe [7] , 
to domain formation in magnetic systems. For an ini- 
tially zero-temperature system, a quench can result from 
rapidly traversing a second order quantum phase tran- 
sition that is associated with a change in the system's 
symmetry. In the case of a quench from SF to MI, the 
SF's U(l) order parameter is absent in the MI phase 
(with its trivial Zi order parameter). In contrast, the 
order parameter transforms from Zi to Z2 for a quenched 
transverse- field Ising ferromagnet [8]. In our experiment, 
we prepare a transversely magnetized two component 
spinor BEC described by a U(l) order parameter, and 
observe the formation and spatial expansion (coarsen- 
ing) of domains following a quench into a phase with a 
U(l) x Z2 order parameter [HHU], unexplored by previ- 
ous studies with binary condensates (miscible [H] [12] or 
immiscible p~3j[14]) with their asymmetric Hamiltonians. 



As compared with three component systems [15-19 , the 
relative simplicity present here allows us to identify an 
intriguing analogy between our spin system and a single- 
component attractive BEC as it collapses [2QH22] . 

Here, we explore the time-evolving magnetization of 
two-component 87 Rb Bose-Einstein condensates (BEC's) 
in the 5Si/2 electronic ground state. Our BEC's 
are well described in terms of a spinor wave-function 
\P(r) = {^(r), ^(r)}, where the |t, I) pseudo- spins label 
|/ = 1,771f = ±1) atomic spin states. The dynamics of this 
wave function are given by the spinor Gross-Pitaevskii 
equation (sGPE) 
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a continuum analog to the transverse field Ising model. 
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In this expression, n(r) 

total density; m is the atomic mass; V(r) is a spin- 
independent external potential; Qj_ describes the Zee- 
man shift of a "transverse" magnetic field; and co,2 
are the spin-independent and spin-dependent interac- 
tion coefficients |23l [24]. In 87 Rb's / = 1 manifold, 
c = (100.86)x47rft 2 a B /m vastly exceeds c 2 ~ -4.7xl0~ 3 c , 
where is the Bohr radius. For a static density profile 
and when Q± = 0, each spin component in Eq. ([!]) is 
separately described by an attractive single-component 
GPE. 

We produce nearly pure N = 7.0(5) x 10 5 atom 87 Rb 
BEC's [28] in the 1/ = 1, = 0) hyperfine state, origi- 
nating from cold |/ = l,ra.p = — 1) thermal clouds formed 
in a hybrid magnetic/optical trap [25]. These BEC's 
are subject to a uniform magnetic field with magnitude 
Bq = 107.0(2) /iT and are confined in the extremely 
anisotropic crossed optical dipole trap [29] depicted in 
Fig. [l^. The radial (e r , i.e., in the e x — e y plane) and 
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FIG. 1: Magnetization M z (r). (a) Schematic representation 
of our experiment, illustrating a spinor BEC with domains 
in an anisotropic crossed-dipole trap, (b) Images showing 
the progression from an initially uniformly magnetized con- 
densate (short times) in which domains appear (intermedi- 
ate times), and then grow spatially (long times); during this 
whole process the condensate slowly decays away, (c) The 
color scale indicates both the degree of magnetization (colors 
from blue to red), and the number density (intensity from 
black to colored). 

axial (e z ) trap frequencies are uj r j2i\ = 135(3) Hz and 
uj z /2n = 3.1(2) Hz respectively. Our T = 90(8) nK con- 
densates have radial and axial Thomas- Fermi radii [30] 
of # r = 3.9(1) fim and = 170(7) fim. 

Because the typical con(r) spin-independent energy 
vastly exceeds the C2n(r) spin-dependent energy scale, 
we make the conventional Thomas-Fermi approximation 
for the overall density distribution n(r) characterized by 
a chemical potential /i, and a minimum healing length £ = 
fi/V2ra/L This gives n(r) = - V(r)] / [c + c 2 M 2 (r)] , 
which depends very weakly on the z component of local 
magnetization vector, M(r) = {M x (r), M y (r), M z (r)} = 

{2Re[V»|(r)^(r)],2Im[^(r)^(r)], |^ t (r)| 2 - |^(r)| 2 }. 

The spin degrees of freedom vary almost exclusively 
with axial position [19] because our extremely anisotropic 
condensate's « 3.9 fim radial extent is comparable to 
the minimum spin healing length £ s = ^"|co/c2 1 1//2 = 
3.20(4) jam. Theoretically, we may describe the spin de- 
gree of freedom as ID spinor with components XT,l( z ) = 
IXt^WI e^,+(*0; retaining terms through first order in 
C2/C0, we obtain an effective ID sGPE 
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The ID interaction strength gir>(z) c-2 is related to a 
ID healing length £ m « y/3/2£ 8 . These two ID sGPE's 



are coupled by the local constraints | ^ (z) | 2 + 1 ^ (z) | 2 = 1 
and (j)^(z) + <f>i(z) =0 (i.e., no mass currents). To make 
the analogy explicit, we have dropped terms quadratic 
in |xt,l| 2 that result from integrating out the transverse 
dimensions. These repulsive terms do not affect the spec- 
trum at short times after the quench, but must be in- 
cluded at long times. 

Our spinor experiment is initiated by a 34 /is rf-pulse 
that puts each atom into a equal-amplitude coherent su- 
perposition of the |t?i) = \ m F — ±1) spin states, the 
ground state when the transverse field Q± is large; the 
system then evolves according to Eq. ([!]) with Q± = 0. 
This procedure is therefore equivalent to rapidly quench- 
ing f2_L to zero: the ground state goes from breaking a 
U(l) symmetry to breaking a different U(l) along with 
a Z2 symmetry (leaving behind a "sneaky" U(l) symme- 
try). The quenched binary mixture is held for a variable 
duration thoid? up to 20 s, while spin structure forms and 
evolves. Spin mixing collisions are suppressed because 
the relatively large 82 Hz quadratic Zeeman shift greatly 
exceeds the C2n(r) ^ 6 Hz spin dependent energy [15]. 
As a result, we observe no population in mp = for 
the entire duration of our experiment. After thoid, we re- 
move the confining potential and allow the atomic ensem- 
ble to expand (largely transversely) for 19.3 ms, during 
which time we Stern-Gerlach [26] separate the spin com- 
ponents. We detect the resulting density distribution by 
absorption imaging, and reconstruct both M x (x,z) and 
M z (x,z), projected onto the e z — e x imaging plane [31]. 
Representative reconstructions of M z (x,y) at six hold 
times are depicted in Fig. ^p. 

The initially (t hold = 0) uniform X (z) = (It) + ID) /V% 
spin superposition is dynamically unstable (as indi- 
cated in Fig. [TJd's snapshots). At this unstable 
point, small spin- wave excitations have an (huj / /lis) 2 = 
(&£id) 2 (&£id) 2 — 2 energy spectrum [5 , where /iiD = 

/i 2 /2ra£ 2 D is a typical ID spin interaction energy. When 
huj is imaginary - for fc£iD £ (0, \/2) - the associated 
modes grow exponentially with peak gain at k = 1/ £id , 
amplifying any existing spin fluctuations, classical or 
quantum. Figure [2] depicts the magnetization M z (z), 
showing the initially unmagnetized condensate develop 
visible structure after about 200 ms. The experimen- 
tal data plotted in Fig. [2^i is in near perfect agreement 
with a a stochastic-pro jective GPE (SP-GPE) simula- 
tion [27] , with parameters nearly matched to our experi- 
ment, Fig. [2}d. 

While the amplitude of these spin waves grow with 
an exponential time constant r(k) = 1/Im(cj(fc, z)) [min- 
imum at t{z) = 2m£ ) 2 B (z)/htt 42 ms], Fig. p] shows that 
no structure is visible until thold ~ 200 ms. Our simula- 
tions confirm that structure begins to grow immediately, 
and only technical noise prevents us from detecting the 
growing spin modulation at shorter times. Figure [2] also 
shows that spin structure forms more slowly in the lower 
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density periphery of the system where £id and therefore 
r are larger. To quantify this effect, Fig. [3] depicts the 
number of spin-regions visible above the noise. This num- 
ber, increasing for short times, is plotted along with the 
results of our SP-GPE simulations, and a local density 
approximation (LDA, accounting for our systems inho- 
mogeneous density profile) prediction for the expected 
pattern of domain growth. 

The spin modulations continue to grow in amplitude 
until, at thoid ~ 300 ms, they form fully spin polarized 
domains of and ||), with a spacing set by the dy- 
namic growth process, not by the system's equilibrium 
thermodynamics. After this period of rapid growth, the 
polarized spin domains evolve slowly, equilibrating, for 
the remaining 20 s duration of our experiment [32] . 

Our BEC has a r = 10(l) s lifetime, implying that the 
domain pattern must evolve in time as the BEC slowly 
contracts. The simplest model - in which the domain 
pattern contracts in lock-step with the dwindling BEC 
- is obviated by Fig. [3J that shows the number of do- 
mains decreasing after thoid ~ 1 s. Indeed, once a domain 
becomes smaller than « 2£id(z), ^ can no longer reach 
full spin-polarization in its center, and it therefore ceases 
to be a barrier for the hydrodynamic flow of the other 
spin state. As a result, small domains de-pin and can 
move freely until they coalesce with another domain of 
the same spin. 

While Figs. |2] and [3] qualitatively suggest that the do- 
mains gradually expand as thoid increases from 300 ms 
to 20 s, it is difficult to obtain a quantitative mea- 
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FIG. 2: Time evolution of magnetization M z (z). (a) Experi- 
mental data and (b) finite temperature simulation using the 
SP-GPE method. In both simulation and experiment, the 
spatial structure of M z (z) coarsens after an initial growth pe- 
riod as domains coalesce. 
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FIG. 3: Number of domains as a function of thoid- The red 
symbols depict the experimentally observed number of do- 
mains (typical uncertainty plotted on the leftmost point) and 
the blue curve plots the results of our SP-GPE simulation 
(uncertainties denoted by the blue band). In both cases, the 
uncertainties reflect the standard deviation over many real- 
izations. In addition, the red curve fits the data to a model 
assuming exponential growth along with a non-zero observa- 
tion threshold, in the LDA. The grey symbols correspond to 
the ratio Rz/^id' an estimate of domain number, assuming 
the system with length 2R Z is partitioned into domains of size 

sure of domain size from data in this form. Indeed, 
the data show that while measurements at neighboring 
times have similar domain sizes, the exact domain pat- 
tern has a significant element of randomness - primar- 
ily in the form of phase shifts - likely resulting from 
subtle differences in the initial conditions, as amplified 
by the subsequent exponential gain process. To miti- 
gate these effects, we turn to the power spectral density 
PSD^ (k) = I J M X:Z (z) ex.p(ikz)dz | obtained from these 
data. With the PSD, we can compare different realiza- 
tions even in the presence of spatial phase shifts of the 
domain structure. 

FigureEk shows PSD z (k) derived from M z (z) shown 
in Fig. [2j For short times (thoid ^ 300 ms), a narrow 
peak associated with the growing spin modulations devel- 
ops. Once the spin domains reach unity polarization, the 
magnetization's magnitude saturates and the boundaries 
between domains - domain walls - sharpen, broadening 
PSD z (k) starting at thoid ~ 300 ms. At longer times, 
the broad peak drifts to smaller wave- vector, indicat- 
ing an increasing typical domain size. Figure com- 
pares this peak location for both experiment and theory 
(red and blue symbols respectively, showing nearly iden- 
tical behavior) against 1/£id- Our simple model predicts 
maximum gain at this wave- vector, however, the peak in 
PSD z (k) for both the experiment and the SP-GPE is at 
slightly smaller k. This results from the inhomogenous 
density profile of the spinor BEC: simulations of uniform 
systems do show peak gain at 1/£id- 

Unlike PSD 2 (/c), PSD^/c) is peaked about zero; this 
is because M x (z) is only appreciable in the domain walls 
where the gas is not fully polarized: it is a series of nar- 
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FIG. 4: Power spectral density, (a) PSD 2 (/c) as a function 
of thoid showing the formation of a peak at finite wave- vector 
k, followed by the gradual movement of that peak to smaller 
k as the spin domains expand. The color scale depicts in- 
creasing spectral power with darker color, (b) Wave- vector 
of PSD z (fc)'s peak, (c) Width of PSD^fc), which always 
peaked around zero. In both (b) and (c), the red symbols 
depict the experimentally observed peak location (typical un- 
certainty plotted on the leftmost point) and the blue curve 
plots the results of our SP-GPE simulation (uncertainties de- 
noted by the blue band). In both cases, the uncertainties re- 
flect the standard deviation over many realizations. The grey 
symbols mark 1/£id, the homogenous-system wave- vector of 
maximum gain (the uncertainties are comparable to the sym- 
bol size). The oscillations for thoid < 1 s in the simulation 
(blue curve, panel c) result from a damped breathing mode 
predominately along e z . 

row peaks. By showing that the width of the peak in 
PSD^fc) tracks the inverse spin-healing length, Fig. [4^ 
demonstrates that the domain walls are sized according 
to fiD (grey symbols). 

Because the ~2£id(z) minimum domain size increases 
as the condensate depletes away, it is plausible that the 
increase in domain-size results from an increasing cut- 
off in the minimum domain size. We exclude this model 
by comparing computed in the BEC's center (Fi 
grey symbols) to the wave- vector of the peak in PSD 2 
Fig. [4] plots experimental data with red symbols and SP- 
GPE simulation with the blue curve. They follow quite 
different time-dependences, with the domains growing 
rapidly at short times before £id appreciably changes, 
indicating that £id is not simply related to the typical 
domain size. 

We conclude instead that we observe genuine coarsen- 
ing of the domains resulting from dynamical exchange of 
particles between stable domains. On the basis of our 
data, we cannot identify the origin of this coarsening, 
but our SP-GPE simulations suggest that spin transport 



through the uncondensed (thermal) fraction, not quan- 
tum tunneling, is the leading mechanism at our temper- 
ature. 

For c 2 < 0, as in 87 Rb, Eqs. and Q allow us to 
describe our system's spin degree of freedom as a single 
component attractive BEC (while the overall density fol- 
lows the conventional Thomas-Fermi profile). In our bi- 
nary mixture, the process of domain formation is a spinor 
analog to the "chain of pearls" pattern that forms in 
ID BEC's quenched from repulsive to attractive interac- 
tions [2Ql[2T] . In that case, the growth of structure results 
from a modulational instability with peak gain at k = l/£ 
set by the conventional healing length. Attractive Bose 
systems are intrinsically unstable against collapse [22 , 
however, in our spinor system, any eventual collapse is 
stymied by an effective hard core interaction resulting 
from the strictly bounded individual spin wavefunctions, 
and higher order interaction terms omitted from Eq. ([2| . 

Thus, in this binary spin system, we observe the full 
gamut of time scales starting with the rapid dynam- 
ical generation of spin-domains from an initially non- 
equilibrium system followed by their subsequent relax- 
ation to progressively larger domains, i.e., coarsening. In 
this model, the coalescence of domains minimizes the at- 
tractive energy by maximizing the spatial extent of re- 
gions with the same spin. Indeed, the ground state con- 
sists of just two domains - one for each spin - thereby 
reducing to one the number of domain walls. The do- 
mains increase in size very slowly in time, but due to the 
overall decrease in the BEC's atom number, we cannot 
distinguish between different functional forms. 
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